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Abstract 

This paper deals with the propagation of Rayleigh waves in 
a homogeneous, prestressed elastic layer of finite thickness 
over a homogeneous, prestressed elastic half-space. The 
dispersion equation has been derived for a layer over a half- 
space, when both media are considered as prestressed and 
the effect of initial stress shown in earlier investigators, is in 
general not applicable to the case of prestressed media, as 
well as the problem of Rayleigh waves has been solved by 
other method without introducing the displacement 
potentials. 

The velocity equations for Rayleigh waves are obtained for 
two cases: 

(I) in a homogeneous, prestressed elastic half-space without 
any superficial layer. 

(II) in a homogeneous, elastic and finite layer over a 
homogeneous elastic half-space. 
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Introduction 

During earthquakes, the Rayleigh waves play more 
drastic role than other seismic waves in damages to 
human beings and buildings at the surface of the Earth. 
Therefore, it is of great importance to the seismologist 
to study the effect of initial stress on the propagation 
of Rayleigh waves. After the pioneering work of 
Rayleigh, many investigators have solved the Rayleigh 
problem for a half-space and one or more superficial 



layers situated over a half-space in homogeneous and 
non-homogeneous media. Most of the works related to 
the Rayleigh waves is well summarized in various 
books i.e. Achenbach, Ewing et al., Miklowitz, Ben- 
Menahem and Singh. However, all these studies have 
ignored the initial stress of the involved media. The 
model of the Earth is supposed to be composed of 
several layer of different thickness and inner layer of 
the Earth under immense stress owing to different 
physical causes i.e. presence of overlying layers, 
variation in temperature and gravitational field. 
Fortunately, Biot developed the incremental 
deformation theory for prestressed medium. Adapting 
the same theory of initial stress, Dey and Mukherjee 
solved the Rayleigh problem for a prestressed layer of 
compressible material and of finite thickness lying 
over an incompressible half-space under hydrostatic 
compressive stresses caused by gravity. The author 
solved this problem by introducing potential method, 
which is, in general, not applicable in prestressed 
media as shown by Sidhu and Singh.The same applies 
to the results obtained by Chattopadhyay et al. for the 
Rayleigh waves problems. In the present paper, our 
aim is to solve the problem firstly by direct method 
then potential method. The frequency equation of 
Rayleigh waves is derived from a homogeneous 
prestressed elastic layer over a homogeneous stressed 
elastic half-space. Here authors deduced two 
particular cases from this main frequency equation. 
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Firstly we obtain a frequency equation for Rayleigh 
waves in a prestressed half-space, which showed that 
this result is similar to the result obtained by Singh 
and Singh for equation of Rayleigh waves 
homogeneous. Secondly we deduced the expression 
for frequency equation of Rayleigh waves in a finite 
homogeneous elastic layer over a homogeneous elastic 
half-space from our main frequency equation by 
putting initial stress equal to zero. The derived result 
is similar to the result given by Ewing et al. 

Basic Equations 

Consider a homogeneous prestressed elastic layer of 
thickness H lying over a homogeneous, prestressed 
elastic half-space. The system is referred to a 
rectangular co-ordinate system with z-axis' directed 
normally downwards and the origin at the interface as 
shown in FIG. 1. Let the principal directions of initial 
stress be chosen to coincide with the directions of 
elastic symmetry and the co-ordinate axes. The state of 
initial stress is, therefore, defined by the principal 
components 5 1X , S 22 and S 33 of the initial stress. We 
restrict our analysis to plane strain parallel to the x — z 
plane with displacements u - win the x- z directions, 
respectively. The third principal stress S 22 = S 33 does 
not enter explicitly into the equations of motion. Its 
influence is, however, included indirectly in the values 
of the incremental elastic coefficients which appear in 
the equations of motion. Let S 1± andS 33 be the initial 
stresses along the x-z directions, respectively in the 
half-space and S^and S 33 be the corresponding 
quantities in the layer. The equations of motion for 
incremental deformations in the half-space are given 
by Biot . 
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Where 



A 1 =Q+ P -, A 2 = Q- P -, 



A 3 =B 12 + Q--, P = S 33 



(2) 



For the layer, we denote the displacement by 

(u,v ,w ) the density by p and the incremental elastic 

coefficients by (fi u , B 22 ,B 12 , Q ). Then we can write the 

equations of motion for the layer similar to equation 

(1). 

Boundary Conditions 
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FIG. 1 GEOMETRY OF THE PROBLEM 

For a layer with free surface at z= -H, the boundary 
conditions are 



Af x =0, A/ z =0atz=-H. 

At the interface between layer and half-space at z = 0, 
the boundary conditions are 

u = u ,w = w , Atz = 

Af x = Af x and A/ z = A/ z . 



(3) 



(4) 



A f x and Af z are the incremental boundary forces per 
unit initial area for the upper layer and the same for 
the half-space being Af x and Af z . The explicit 
expressions for Af x ,Af z ,Af x ea\d Af z are 



Af x — S 31 + S 33 a> y S\\£zx> 
A/z = + S 33 e xx , 



(5) 
(6) 
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and 



(7) 



4/z — ^33 + ^33 e %% -$31 > ^33 '^31 '^33 are 

incremental stresses, e zx , e xx , e zx , e xx are incremental 
strainsand a> y , a> y are the incremental rotation 
components parallel to the xz-plane. Explicit 
expressions of these quantities in terms of u and w 
are 



du , „' dw 



S 33 - (#12 - P ) -T7 + B 2 



dx ■ ~ 22 dz ' 



(8) 
(9) 



, _3u' 
6xx ~ ~dx' 

1 /au' aw'^ 



(10) 



(ii) 



e ' zx ~2\dx + dz J' 

"y-iV^-^z-)' 

with similar expression for S 33 etc. 
Propagation of Surface Waves 



For Rayleigh waves of circular frequency a>, wave 
number k and phase velocity C R propagating in the x- 
direction through the half-space, we may assume the 
solution of equation(l) 



u = u r e~ kqz e ik ( CRt ~ x ^ 
w = w r e~ kqz e ik ^ Rt ~ x ^ , 



(12) 



where U r and W r are the amplitude factors and q is 
assumed to be real and positive. Inserting the values of 
the displacements in equation (1), we get 



(pC| - B u + A ± q 2 )U R + iqA 3 W r = , 



(14) 



= 0. 



iqA 3 U r + (pC 2 +B 22 q 2 -A 2 )W r 
(15) 



The values of q may be obtained from 
18 



(pC 2 - fi n + A iq 2 ) iqA 3 

iqA 3 (j)Ci + B 22 q 2 - A 2 ) 

= 



On which expansion takes the form 

? ? 

- [A 2 + pC 2 ^ + B 22 ) - B 1± B 22 - A ± A 2 ] + 
[{A 2 + pC 2 ^ + B 22 ) - B 1± B 22 - A ± A 2 } 2 - 
= 4{A 1 B 22 (A 2 B 11 + p 2 Ct - pC 2 A 2 - pCjB^}] 1 ^ 
2A 1 B 22 

therefore the solution (13) can be written as 

u = (U ri e-^ kz + U r2 e-^ kz )e ik ^ CRt - x \ 



(16) 



(17) 



w 



= (W ri e~^ kz 

+ W r e -l2kz^ e ik(C R t-x) _ 



(18) 



(19) 



In the solution (13) of the equation of motion (1) for 
the half-space, we have considered only e~ qikz so that 
the radiation condition is satisfied. However, for the 
layer we must take both e~ kq z and e kqz terms. For 
e -kq z tne so i u ti on in layer is similar to equation (19). 
However, for e kq z , we may assume 



= U t e kq z e ik ( c R t - x \ 



w 



(13) = w t e kq z e lk ^ R t-x)_ 



(20) 



(21) 



Inserting the values of the displacements in equation 
(1), we get 



(p'c 2 - fi n + A iq 2 ) U't ~ iq'A 3 W t = 0, 



(22) 



-iq'A 3 U t +(pC 2 +q B 22 - A 2 )W t = 0. (23) 
The values of q may be obtained from 
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p'CZ-B[ 1+ AW 



-iq'A' 3 



-iq'A' 3 
p'C 2 +q' 2 B 22 -A' 2 



= 



(24) 



and 



q A1B22 + <rUi + pcUi + pc£B 22 - b x1 b 22 - 

A ± A 2 ) + (A 2 B 1± + p 2 C*-pC£A 2 - pC|Bn) = . (25) 
This is a quadratic equation in q 2 and whose roots are 

,2 ,2 9 9 

q 1 , q 2 similar to the values of q{ and q 2 as in 
equation (17). Therefore the solution (19) can be 
written as 



= {U' ti e k ^ z 

+ U' t e k i'i z ^ e ik VR t - x \ 



(26) 



w 



Surface for the upper layer 



(27) 



by equation (15) for q = q 1 and q 2 . Taking second 
member of equations (23) and (15), we get 

(p'C| + q'l B 22 - A 2 ) W tl = , (30) 

-iq 2 A' 3 Ul 2 + 

{p'Cl+q' 2 B' 22 -A' 2 )Wl 2 = Q, (31) 
iq[A' 3 U ri + {p'Cl + q[ 2 B' 22 - A' 2 )W r [ = , (32) 
iq' 2 A' 3 U; 2 + (p'Cj + q' 2 B' 22 - A' 2 )W r ' 2 = , (33) 



Equation (33) may be written as 
U' r ^=m\W^ ,u; 2 =-m' 2 W r ' , 

U't^m!^ ,Ul 2 = -m' 2 w; 2 , 
where 

= iM[,m' 2 = iM' 2 , 

and 



(34) 



(35) 



(36) 



We may, therefore, assume the total displacement field 
to be the form in the layer 



W = 



(U' t± e k ^ z + U' t2 e k ^ z + U r± e~ k ^ z + 

\]' r e -'"?2Z) e ifc(CRt-3C) j 



(28) 



w 



= [W^e kq '^ z + Wl 2 e kq '^ z + W^e~ k ^ z 

+ W r ' 2 e- kq z z) e ifc(c R t-^) j (29) 

uQ , (w? lf wQ , {u' ri , u; 2 ) , (w; it w; 2 ), 

C R and k have conventional 

meaning in the layer. q 1 , q 2 are real and positive 
quantities. 

(U tl , U t2 ) and (W tl , VK t2 )are not independent but 
connected by equation (23) for q = q 1 and q 2 . 
Similarly {u' Tl , U r2 ) and (W r± , W ri ) are also connected 



p'C 2 +q[ 2 B 22 -A' 2 

Ml = 541 



, _ p'Cj + q' 2 B' 22 - A' 2 



M[ = 



q' 2 A' 3 



(37) 



Solution for the half-layer 

The displacement field in the half-space may be 
written in the form equation (11). 



= {U'r 



-kq t z 



= {Wr x e~ k ^ z 

+ Wy e ~ k( l2^ e ik{C R t-x) _ 



(38) 



w 



(39) 
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(t/ ri , V r2 ) and (W ri , W r2 ) are not independent but 
connected by equation (15) for q = q 1 and q 2 . Taking 
second member of equation (15), we get 

i qi A 3 U ri + {pC 2 R + qlB 22 - A 2 )W ri = , (40) 

iq 2 A 3 U r2 + (pC 2 R + q\B 22 - A 2 )W r2 = . (41) 
Equation (41) may be written as 

U ri = mi W ri ,U r2 =m 2 W r2 , (42) 
where 

m 1 = iM 1 , m 2 = iM 1 (43) 
and 



M, = 



PCr +qlB 22 -A 2 



qiA 3 



My = 



pCj + q\B 22 - A 2 



(44) 



(U ri , \] r2 , W ri , W r2 ), C R and k have conventional 



meaning in the half-space. q 1 , q 2 are real and positive 
quantities defined in equation (17). The total 
displacement field given in equations (29) and (39) 
must satisfy the boundary conditions equations (3) 
and (4). Making use of equation (34) and (42), we 
obtain 



ae- k ^"W tl + be~ k ^ H Wi 2 + ae~ k ^ H Wr x + 



b e - k ^ H Wr 2 + 0W ri + 0W ri = , 



(45) 



f e - k ^"W tl + ge- k ^"W t2 - fe~ k ^ H W ri - 

ge - k ^ H W r2 + 0W ri + 0W r2 = 0, (46) 

aW^ + bW{ 2 + aW; t + bw; 2 + lW ri + mW r2 = , (47) 

fW^ + gWl 2 - fW r [ - gW; 2 + nW ri + rW r2 = , (48) 



-m^W^-m^W^ + m' x w; x + m! 2 W^ - mi W ri - m 2 W r2 = 
, (49) 
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Wt t + K + w; t + w; 2 - w ri -w r2 = o, 

where 



P'\ ( R' 

a = - I Q '+y Jmiqi + I Q ~~2 I' 



(50) 



m 



Q'+Y^rn 2 q 2 + (Q'-^-)i 

f = {B' 12 +S' 11 )m' 1 i + B 22 q[ , 
g = (B[ 2 + S^m^i + B' 22 q 2 , 

= [{Q+^m iqi + (Q-^i 

= [(Q+^)m 2q2 + (Q-^i 
f = (B 12 + S^m^ + B 22 q u 



(51) 



r = (B 12 l-S^m^ + B 22 q ± 

Eliminating W tl , W t2 , W ry W r2 , W Tl and W rr we obtain 
the following frequency equation for Rayleigh waves 
in prestressed media. 

DetA 

= (52) 
where A 
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Particular Cases 

If the initial state is unstressed, the elastic coefficients 
for the half-space will become 

fin = B 22 = X + 2fi, A 1 =A 2 = Q= t i, 
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B 12 =A, A 3 =A + n,R = P = 0. 
From equations (44) and (17), we get 



M 1 = q 1 ,M 2 = — , 

R2 



(53) 

(54) 
(55) 



respectively, where 



, A + 2/i 
cr = ana p £ 



P ' 



A and ji are Lame's constant for the half-space in 
unstressed media .Similarly for the layer, 
corresponding quantities becomes 

fin = #22 = A' + 2// ,4 = A 2 = Q' = n' 

B 12 =A,A 3 = A + [i ,R = P = , 
and 



(56) 



(57) 



M 1 = q v M 2 = — , 

<?2 

where 



(58) 



(59) 



Using equations (53) - (59) in equation (51) becomes 



a = — i I 2 7 M -b = —2iu , 



.1 "J J /' ■.'/ 



V 2 ' 



l = i\2- — )n ,m = 2ifi, 



n = 2(l-4) '\, r = 



(60) 



Rayleigh Waves in Prestressed Half-space 

In the absence of any prestressed layer, only two 
boundary conditions will be satisfied by the solution 
of the half-space. These are Af x =Af z = 0. Eliminating 
W ri and W r2 , the frequency equation may be written as 
1 1 ml 
In r I 



0, 



Which after simplification of whichthe form can be 
obtained 

[(<? + ^)rn 1 qi + i ( Q - [i(B 12 + S 11 )m 2 + B 22 q 2 ] - 



[i(B 12 + S^nh + B 22qi ] [(Q + m 2 q 2 +l(Q- |)] 

(61) 



= 



This result is equivalent to the corresponding result of 
Singh and Singh for the frequency equation of 
Rayleigh waves in the prestressed, homogeneous 
elastic half -space. 

Rayleigh Wave in Unstressed Media 

when both media are unstressed i.e. P = R = P = R = 
0. Making use of equations( 53)-(36), the equation (52) 
takes the form 



DetB = 

where B is 



(62) 



which is the frequency equation of Rayleigh waves in 
a layer over a half-space in the unstressed media. 
These results coincide with the results given by Ewing 
et al.. 

Conclusions 

For a prestressed layer lying over a prestressed half- 
space, the propagation of Rayleigh waves is 
considered. It is noted that the frequency equation of 
Rayleigh waves are effected due to the initial stresses 
present in the equation. The particular result is 
compared with the standard result 
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Det B 



(-1 



a' I 



2(1 



a' 



2-% 



V 2 



r 2 



(1--t) 



a' 



r 2 



r 2 
a' 



2 -«).-H*) v 



r 2 



-(l--r) 



f Cr 
— (2--^) 
/?' 2 



(l-^)V 2 



2 £(2-4)V2 



obtained by Singh and Singh for frequency equation in 
prestress half-space. The frequency equation of 
Rayleigh waves for homogeneous elastic layer over a 
homogeneous elastic half-space is obtained and 
compared with the result without stresses given by 
Ewing et al.. 
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